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The body temperature is considered a universal cue by which the master clock synchronizes the peripheral 
clocks in mammals, but the mechanism is not fully understood. Here we identified two cold-induced 
RNA-binding proteins (RBPs), Cirbp and Rbm3, as important regulators for the temperature entrained 
circadian gene expression. The depletion of Cirbp or Rbm3 significantly reduced the amplitudes of core 
circadian genes. PAR-CLIP analyses showed that the 3'UTR binding sites of Cirbp and Rbm3 were 
significantly enriched near the polyadenylation sites (PASs). Furthermore, the depletion of Cirbp or Rbm3 
shortened 3'UTR, whereas low temperature (upregulating Cirbp and Rbm3) lengthened 3'UTR. 
Remarkably, we found that they repressed the usage of proximal PASs by binding to the common 3'UTR, 
and many cases of proximal/distal PAS selection regulated by them showed strong circadian oscillations. 
Our results suggested that Cirbp and Rbm3 regulated the circadian gene expression by controlling 
alternative polyadenylation (APA). 
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The circadian clock is a daily time-keeping mechanism that regulates a wide range of biochemical, physio- 
logical and behavioral processes in living organisms'. In mammals, most cells in the body contain cell- 
autonomous circadian clocks 2,3 . The circadian clocks of peripheral tissues were synchronized by the master 
clock residing in the suprachiasmatic nucleus (SCN) 4 . In this synchronization process, the circadian oscillation of 
body temperature controlled by the SCN, can serve as a global entrainment cue 5,6 . The oscillation of body 
temperature consists of the temperature-rising phase and the temperature-declining phase, which can be con- 
sidered a heat shock process and a cold shock process, respectively. The heat shock pathway was proposed to play 
a critical role in synchronizing the peripheral clocks to thermal stimuli 5,7 . However, the modulators in the 
temperature-declining phase remain to be identified. Here we proposed that cold shock proteins, such as cold- 
induced RBPs, may regulate the circadian clock in this process. 

RBPs have long been known to play important roles in the cellular response to low temperature. In E. coli, the 
CspA can destabilize the RNA secondary structures during hypothermia and thereby facilitate translation 8 . The 
Y-box binding proteins in Xenopus oocytes that contain a cold-shock domain can control translation during 
hypothermia 9 . In mammals, two cold-induced RBPs, Cirbp and Rbm3, were shown to function as RNA chaper- 
ones and influence both transcription and translation during hypothermia 1012 . Cirbp and Rbm3 belong to the 
same family of cold- induced RBPs and are highly homologous to each other 10 . Although it was initially discovered 
in response to the UV irradiation 13 , Cirbp was later found to be induced by other stressors as well, such as 
hypothermia and hypoxia 14,15 . Rbm3 is one of the most significantly upregulated genes in multiple tissues during 
animal hibernation 16 . 

Increasing studies indicated the importance of posttranscriptional regulation by RBPs in the mammalian 
circadian clock, but the regulatory mechanisms have been limited in the mRNA translation and stability control 17 . 
Rbm4 (mLark) can promote the translation of Perl by binding to the 3'UTR 18 . The RBPs Ptb and Hnrnpd were 
shown to promote the degradation of Per2 and Cryl mRNAs by binding to their 3'UTRs, respectively 19,20 . In 
addition, most recently, Cirbp was shown to modulate circadian gene expression posttranscriptionally, but the 
detailed mechanisms were not fully elucidated 21 . In this study, we proposed a novel posttranscriptional regulatory 
mechanism in circadian clock: APA regulation. 



SCIENTIFIC REPORTS | 3 : 2054 | DOI: 1 0. 1 038/srep02054 



1 



APA is one of the most important mechanisms in the posttran- 
scriptional regulation. More than half of the mammalian genes con- 
tain tandem UTRs which are generated by multiple PASs in the 
terminal exon 22 . Without changing the mRNA coding potential, 
tandem UTRs can alter the length of 3'UTR, thus influence the 
mRNA stability, translation and localization by changing the access- 
ibility of RBPs and miRNAs to 3'UTR 23 . Previous studies have shown 
that 3' UTRs tend to be shortened during T lymphocytes activation 24 
and in cancer 25 , but lengthened during embryonic development 26 . 
However, the regulators of APA during specific biological processes 
remain to be identified. 

In this study, we demonstrated a new mechanism in the temper- 
ature entrainment of circadian clock. We showed that two cold- 
induced RBPs (Cirbp and Rbm3) can control the APA of their target 
genes and influence the amplitude of temperature entrained cir- 
cadian clock. Cirbp and Rbm3 were required for high- amplitude 
oscillation of core circadian genes in the temperature synchronized 
cells. The 3'UTR binding sites of Cirbp and Rbm3 were enriched 
within 100 nucleotides upstream of the PASs, and these two RBPs 
repressed the usage of proximal PASs by binding to the common 
3'UTR. Furthermore, we found that the usage of proximal or distal 
PASs which was regulated by Cirbp or Rbm3 showed strong cir- 
cadian oscillations in mouse liver. These results suggest that the 
regulation of APA by Cirbp and Rbm3 is an important process in 
regulating the amplitude of temperature entrained circadian clock. 



Results 

Cirbp and Rbm3 are essential for the temperature entrained circa- 
dian gene expression. To identify the candidate genes which are 
crucial in the temperature entrainment of circadian clock, we 
searched for the genes that showed circadian oscillating expression 
in multiple tissues, and their oscillating expression should be depen- 
dent on body temperature and independent of the local molecular 
clock. In a meta-analysis of circadian microarray datasets covering 14 
mouse tissues, 9,887 known genes exhibited circadian oscillations in 
at least one tissue, and 1 5 1 of them showed circadian oscillations in at 
least six tissues 27 . We divided the 151 genes into six groups based on 
their circadian peak times (Fig. la) and performed gene ontology 
(GO) enrichment analysis for each group using DAVID 28 (Table 1). 
The well-known core circadian genes were significantly enriched in 
the CT10 and CT22 phases. Besides, the protein-folding related 
genes encoding heat shock proteins (HSPs) were significantly 
enriched in the CT18 phase (i.e. the temperature rising phase). 
The genes of RBPs were significantly enriched in the CT6 phase 
(i.e. the temperature declining phase). This result was consistent 
with the pre-existing knowledge of HSPs and RBPs functioning in 
response to high and low temperature, respectively 7,14,29 . 

We then analyzed a published circadian microarray dataset of 
mouse liver (ArrayExpress, E-MEXP-842), in which the liver 
molecular clock was either active or arrested 30 . This analysis iden- 
tified 401 genes that exhibited a highly similar circadian oscillating 
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Figure 1 | RBPs are potential regulators in the temperature entrainment of circadian clock, (a) Histogram of the circadian peak times of 151 circadian 
oscillating genes. Black curve represented the mouse daily body temperature 52 . The groups of the genes with significantly enriched GO categories were 
highlighted. CT, circadian time in hour (h). (b) The schematic diagram of the square wave temperature synchronization model. Sampling time was 
indicated, (c-e) The expression patterns of the core circadian genes and the genes of HSPs and RBPs in the temperature synchronization model. The 
lowest expression value was normalized to 1 (black, NN condition; red, NC condition). Data are represented as mean +/— SEM (n = 3). 



SCIENTIFIC REPORTS | 3 : 2054 | DOI: 1 0. 1 038/srep02054 



2 



Table 1 | 


Significantly enriched GO terms 


in each group of genes in Figure 1 A 




rnase 


OO terms 


Oenes 


Fisher s exact test P 


CT6 


RNA binding 


Cirbp, Fus, Hnrpdl, Rbm3, Rbmsl, Srsf5, Srsf6 


1.3E-04 


CT10 


rhythmic process 


Dec), Dbp, Perl, Pert, lex, Nrld2 


2.7E-04 


CT18 


protein folding 


Hspa5, Hspa8, Hspalb, Hsphl, Hspbl, Hsp90aal, Hsp90abl 


1 .9E-06 


CT22 


rhythmic process 


Bmall, Clock, Nfil3, Bcl2 


3.8E-04 



expression between the liver clock active and the liver clock arrested 
mice (Supplementary Fig. Sla; COSOPT 31 , P < 0.01; Pearson cor- 
relation > 0.5). For these genes that showed independent oscillating 
expression of local molecular clock, they must be driven by systemic 
signals, such as temperature. Almost all of the RBPs and HSPs in 
Table 1 belonged to this gene set, which showed persistent oscillating 
expression in the liver clock arrested mice. (Figs. Sib and Sic). This 
result implied that the RBPs and HSPs are potential candidates to 
play important roles in the temperature entrainment of peripheral 
clocks. Previous studies have identified HSPs as important modula- 
tors in resetting the peripheral clocks to thermal stimuli 5 , but the 
function of RBPs in this process needs further investigation. 

To identify the regulator in the temperature entrainment of cir- 
cadian clock, we subjected mouse embryonic fibroblasts (MEFs) to a 
square wave temperature cycle between 37°C and 33°C, i.e. 12 h at a 
normal temperature (N) of 37°C followed by 12 h at a low temper- 
ature (C) of 33°C, to model the temperature entrainment naturally 
occurs in mice 32 (Fig. lb). After six days of temperature cycling, we 
either changed the condition to constant temperature of 37°C (NN 
condition) or continued the temperature cycling (NC condition). 
The expression of the core circadian genes and the genes of HSPs 
and RBPs were measured by quantitative real-time PCR (qPCR) 
every four hours for two days. The expression patterns of the can- 
didate genes can be divided into three groups. The core circadian 
genes such as Bmall and Perl oscillated in both NN and NC condi- 
tions (Fig. lc). Once entrained by temperature cycling, their express- 
ion oscillated no matter the temperature cycling still exists or not. 
The genes encoding HSPs did not oscillate in either NN or NC 
condition (Fig. Id). This result supported our earlier claim that the 
HSPs responded to high temperature but not low temperature. The 
two genes encoding cold-induced RBPs, Cirbp and Rbm3, oscillated 
only in the NC condition but not in the NN condition (Fig. le). Their 
expression levels increased following the decrease in tempera- 
ture and stabilized as soon as the temperature became constant. 
Because the expression of Cirbp and Rbm3 were highly correlated 
with low temperature, oscillating in multiple tissues and independent 
of the local molecular clock, we proposed that the two RBPs played 
an important role in the temperature entrainment of circadian clock. 

Our previous microarray analysis reported that the mRNA levels 
of Cirbp and Rbm3 oscillate in multiple mouse tissues 27 . We exam- 
ined their mRNA levels in mouse liver and cortex (See Experimental 
Procedures), and found that their expression showed strong circadian 
oscillations in both tissues (Supplementary Fig. S2a). Moreover, the 
protein levels of Cirbp and Rbm3 in the liver showed strong oscil- 
lation with peaks around the CT6 and CT10, when the mouse body 
temperature was low (Supplementary Figs. S2b and S2c). 

To examine whether Cirbp and Rbm3 are required for the tem- 
perature entrained circadian gene expression, we employed two 
independent siRNAs to deplete Cirbp or Rbm3 in MEFs (Figs. 2a 
and 2b). We entrained mock, control, Cirbp- and Rbm3-depleted 
cells with the NN condition of temperature synchronization model. 
We monitored the expression of Dbp and Per2 at a series of circadian 
time points. The oscillating amplitudes in the mRNA levels for both 
genes were significantly reduced in both Cirbp- and Rbm3-depleted 
cells comparing with the mock and control cells (Fig. 2c and 
Supplementary Fig. S2d). To determine the regulatory extent of 
Cirbp and Rbm3 on circadian gene expression, we examined the 



expression of core circadian genes at two time points (16 h and 
28 h after temperature synchronization) in the temperature synchro- 
nized cells. Four out of 12 core circadian genes, including Perl, Dbp, 
Nrldl and Nrld2, showed significantly reduced expression changes 
between the two time points in both Cirbp- and Rbm3-depleted cells 
comparing with the mock and control cells (Fig. 2d and 
Supplementary Fig. S2e). These results indicated the importance of 
Cirbp and Rbm3 for regulating the amplitude of the temperature 
entrained circadian gene expression. 

Identification and analyses of Cirbp- and Rbm3-binding sites 
using PAR-CLIP. To study the functions of Cirbp and Rbm3, we 
used PAR- CLIP (photoactivatable ribonucleoside enhanced cross- 
linking and immunoprecipitation) 33 35 in combination with next- 
generation sequencing to identify their direct binding targets. In 
PAR-CLIP, the photoactivatable nucleotide (4-thiouridine, 4SU) is 
incorporated into nascent RNA and thereby enhances the 
crosslinking of protein to RNA. Moreover, the crosslinking of 
4SU-labeled RNA to proteins leads to specific T-to-C changes in 
cDNA sequences, marking the protein binding sites on the target 
RNA 33 , which facilitates the discrimination of the signal-to-noise 
ratio. The Cirbp- or Rbm3-bound RNA fragments were success- 
fully recovered (Supplementary Fig. S3a), purified and converted 
to cDNAs after the ligation of adapters and then sequenced on 
Illumina GAII platform. 

We observed diagnostic T-to-C transitions in PAR-CLIP cDNA 
libraries of Cirbp and Rbm3 (Supplementary Figs. S3b and S3c). 
These T-to-C changes indicated the existence of the bona fide RBP 
binding sites on the target RNA 33 . All of the sequencing data were 
analyzed with a recently developed computational pipeline 34 . Two 
independent PAR-CLIP experiments yielded 23.6 million and 4.98 
million reads for Cirbp, 13.1 million and 12.6 million reads for Rbm3, 
which were aligned to the mouse mm9 genome by allowing only 
unique genomic hits and one nucleotide mismatch (Supplementary 
Table SI). Using the aligned reads, we generated 12,070 and 16,893 
Cirbp- and Rbm3-binding clusters with T-to-C transitions and 20 or 
more overlapping reads, respectively. The mean lengths of the 
Cirbp- and Rbm3-binding clusters were 25.9 and 23.1 nucleotides, 
respectively, which indicated the high resolution of PAR-CLIP. 8,597 
Cirbp- and 9,669 Rbm3-binding sites can be assigned to annotated 
gene regions, including 5'UTR, coding sequence (CDS), 3'UTR and 
intron, based on the UCSC mouse gene annotation 36 (Figs. 3a and 3b, 
Supplementary Table S2). The majority of the Cirbp- or Rbm3-bind- 
ing sites were located in the exonic regions, especially the CDS and 
3 ' UTR. Furthermore, we found that the 3 ' UTR binding sites of Cirbp 
and Rbm3 were significantly enriched within 100 nucleotides up- 
stream of PASs (Cirbp, 939 out of 3,520; Proportion test, P < 1.0E- 
15; Rbm3, 549 out of 2,081; P < 1.0E-15; Fig. 3c), which implied the 
potential function of the two RBPs in the polyadenylation. 

We employed the MEME software 37 to identify the sequence 
motifs in the Cirbp- or Rbm3-binding clusters which were located 
in the exonic regions. In the 7,139 Cirbp-binding clusters, MEME 
identified a strong motif that was present in 3,353 (47.0%) clusters 
(Fig. 3d). In the 6,106 Rbm3-binding clusters, a strong motif was 
present in 3,533 (57.9%) clusters (Fig. 3e). These results indicated 
that the binding motifs of Cirbp and Rbm3 are not identical, but very 
similar to each other. Previous studies have indicated that functional 



SCIENTIFIC REPORTS | 3 : 2054 | DOI: 1 0. 1 038/srep02054 



3 



I Cirbp t Rbm3 

I 
** ** o-H ■ 
ll o I -- 



a 1 ' 

_ o 

<D < 

| |o. 5 - 

E 2 0 j ■ o- 

control si-1 si-2 control si-1 si-2 



Actb 
Cirbp 

Actb 
Rbm3| 



control si-1 si-2 



control si-1 si-2 



30- 



■D CD 



CO < 15- 

11 



Dbp Per2 



i — i — i — i — i — i — i — i — i — i — i — i 

4 12 20 28 36 44 
Time after temperature 
synchronization (h) 



i — i — i — i — i — i — i — i — i — i — r~ 
4 12 20 28 36 44 
Time after temperature 
synchronization (h) 



0 

CD 



o _ 
o 



Perl 



Nr1d1 



CO CN^ J 
CO — 
P r- 



g-CD 
LU C 



0 



12 



Nr1d2 



* * 6- 

lOo. 




mock 
control 
siCirbp-1 
siRbm3-1 



Figure 2 | The depletion of Cirbp or Rbm3 reduced the oscillating amplitude of core circadian genes in the temperature synchronized cells, (a and b) 

The mRNA and protein levels for the control and the two independent siRNA knockdowns of Cirbp and Rbm3. Actb was used as control, (c) The 
oscillating mRNA levels of Dbp and Per2 in the temperature synchronized mock (black), control (gray), siCirbp-1 (blue) and siRbm3-l (red) cells. The 
lowest expression value was normalized to 1. See also Fig. S2C. (d) The expression changes between 16 h and 28 h after temperature synchronization of 
the core circadian genes (Perl, Dbp, Nrldl and Nrld2) in the temperature synchronized mock, control, siCirbp-1 and siRbm3-l cells. The statistical 
significances were calculated between the control and Cirbp- or Rbm3-depleted cells. See also Fig. S2D. Data are represented as mean +/— SEM (n = 3); 
Student's f-test * P < 0.05, ** P < 0.01, *** P < 0.001. 



RNA binding elements tend to be more conserved than their flanking 
genomic regions 34,38 . We therefore estimated the evolutionary con- 
servation of Cirbp- or Rbm3-binding sites across 20 placental mam- 
malian species using the PhyloP nucleotide conservation scores 39 . In 
each binding cluster, the position with the highest frequency of T-to- 
C transitions was considered the preferred crosslinking site. We 
extracted the PhyloP scores at a given distance around the cross- 
linking site and averaged these over all binding sites (Fig. 3f), and 
found that the sequences around the crosslinking site showed highly 
elevated conservation, which suggests the existence of functional 
elements. 

In total, we identified 4,376 and 5,300 target genes of Cirbp and 
Rbm3, respectively. Despite their different binding motifs, their tar- 
get genes had significant overlap (2,539 genes; Fisher's exact test, P < 
1.0E-15). GO enrichment analyses showed that their target genes 
were mostly involved in regulating the transcription and mRNA 
metabolic processes (Supplementary Table S3). In addition, we inde- 
pendently performed Cirbp or Rbm3 RNA-IP in MEFs and validated 
five Cirbp- and Rbm3-target genes out of five tested through qPCR, 
including the core circadian genes Nrld2 and Tef (Supplementary 
Fig. S3d). 

To evaluate the function of Cirbp and Rbm3 in the circadian clock, 
we compared their target genes with the circadian gene database 27 
and found that the circadian oscillating genes are significantly 
enriched. 2,458 (56.4%; Chi-square test, P < 1.0E-15) Cirbp-target 
genes and 2,762 (52.1%; Chi-square test, P < 1.0E-15) Rbm3-target 
genes exhibited oscillating expression in at least one tissue. When we 
compared oscillating Cirbp- or Rbm3-target gene counts to the 



expected oscillating gene counts by chance, we found the enrichment 
ratios were 2.1 and 2.2 for the genes oscillating in one tissue, and the 
enrichment ratios increased with the number of tissues in which the 
genes oscillated (Fig. 3g). We also found that the circadian peak times 
of the Cirbp- or Rbm3-target genes were significantly enriched in the 
CT6 phase (i.e. the temperature declining phase; Supplementary 
Figs. S3e and S3f). Furthermore, many core circadian genes, includ- 
ing Bmall , Clock, Perl , Cryl , Nrl d2 and Tef, were shown to be bound 
by Cirbp or Rbm3. The Cirbp- or Rbm3-target genes that showed 
oscillating expression in at least eight tissues were shown in 
Supplementary Table S4. These results suggested that the two 
RBPs are important regulators in the circadian clock. 

Cirbp and Rbm3 repress the usage of proximal P ASs by binding to 
the common 3'UTR. To assess the functions of Cirbp and Rbm3 at 
transcriptome level, we applied RNA-seq to mock cells, Cirbp- 
depleted cells, Rbm3 -depleted cells and mock cells that were 
treated with cold shock (32°C) (Supplementary Table S5). The 
analyses of the RNA-seq datasets showed that the results between 
the biological replicates were highly reproducible (Supplementary 
Figs. S4a-d). Using stringent criteria (FPKM > 1 in every sample, 
and log2 [fold change] > 0.8), we identified 876 upregulated and 454 
downregulated genes upon Cirbp depletion, and 237 upregulated 
and 293 downregulated genes upon Rbm3 depletion (Supplemen- 
tary Table S6). The GO enrichment analyses showed that the genes 
which were upregulated upon Cirbp or Rbm3 depletion were 
significantly enriched in cell cycle (Supplementary Table S7). This 
result was consistent with the previous studies that Cirbp and Rbm3 



SCIENTIFIC REPORTS | 3 : 2054 | DOI: 1 0. 1 038/srep02054 



4 



Cirbp 



5'UTR, 1% 



\ 



Intron, 
17% 

3'UTR, 
41% 



CDS, 
41% 



Rbm3 



5'UTR, 2\ 

Intron, 
37% 



CDS, 
39% 



3'UTR, 
22% 



160 



£ 
c 
=> 
o 
0 

L_ 

B 

to 

o 



80 



— Cirbp 

— Rbm3 




-500 -300 -100 0 
Distance relative to PAS [nt] 



Cirbp 



1.1E-316 



(l (l 



1 2 3 4 5 6 7 8 



2- 
0. 



Rbm3 2.0E-100 



G 



1 2 3 4 5 6 7 



1.2 



a> 

o 
o 
in 

q_ 0.9 
_g 

>, 



0.6 



Cirbp 



control for Cirbp 



control for Rbm3 



-100 0 100 

Distance from crosslink site [nt] 



6l 

o 
to 

~ 41 
c 
a> 
E 

5 2 



Cirbp 
Rbm3 



* * * * * * 

* * * * * * 

* * * * * * 



* * * * 

* * * * 

* * * * 



* * * * 

* * * * 

* * * * ■-- 




1 2 3 4 5 6 >=7 
Tissue counts 



Figure 3 | Identification and analyses of Cirbp- and Rbm3-binding sites using PAR-CLIP, (a and b) The distribution of Cirbp- and Rbm3-binding sites 
in gene region revealed predominant binding within CDS and 3'UTR. (c) The distribution of Cirbp- and Rbm3-binding sites relative to thePASs revealed 
the enrichment within 100 nucleotides upstream of the PASs. (d and e) The enriched sequence motifs of Cirbp- and Rbm3-binding clusters which are 
located in the exonic regions were identified by the MEME software. For each motif we indicated the E-value. (f) The conservation profiles of Cirbp- and 
Rbm3 -binding sites that are located in the exonic regions. Average PhyloP score at a given distance from the crosslinking sites was shown. Random 
selected regions from the same transcript were drawn to serve as control, (g) Enrichment ratios of Cirbp- or Rbm3-target circadian oscillating genes 
compared to the expected circadian oscillating genes by chance in the certain tissue counts. Chi-square test ** P < 0.01, *** P < 0.001. 



regulated cellular proliferation 40 " 43 . To assess the influence of Cirbp 
or Rbm3 on the stability of their target transcripts, we compared the 
mRNA levels from knockdown and mock cells. The expression levels 
of transcripts with Cirbp-binding sites in the intron significantly 
decreased upon Cirbp depletion, but the levels of transcripts with 
Cirbp-binding sites in the CDS and 3'UTR were not significantly 
changed upon Cirbp depletion (Supplementary Fig. S4e). The 
expression levels of transcripts with Rbm3 binding sites in CDS, 
3'UTR and intron significantly decreased upon Rbm3 depletion 
(Supplementary Fig. S4f). These results suggested that Cirbp 
stabilized its target pre-mRNAs, and Rbm3 stabilized its target pre- 
mRNAs and mRNAs. 

Because the PAR-CLIP results implied the potential functions of 
Cirbp and Rbm3 in the polyadenylation, we systematically screened 
our RNA-seq data to identify the role of Cirbp and Rbm3 in the APA 
regulation. After combining the PolyA_DB2 database 44 with the lat- 
est polyA-seq datasets 22 , we annotated 12,250 genes with reliable 
tandem UTRs from the UCSC known genes in mouse. We then 
compared the relative expression of the common 3'UTR (between 
stop codon and proximal PAS) and of the extended 3'UTR (between 
proximal PAS and distal PAS) upon Cirbp or Rbm3 depletion, and 
identified the switching events in the usage of proximal or distal PAS 
(PAS-usage switching events) which were regulated by Cirbp or 
Rbm3. We identified 247 and 140 genes that showed significantly 
changed ratios of extended to common 3'UTR (ext/com UTR-ratio) 
upon Cirbp or Rbm3 depletion under the criteria (Fisher's exact test, 
P < 0.01; the ext/com UTR-ratio change > 0.6). Moreover, 166 



(67.2%) and 119 (85.0%) of the genes showed downregulated ext/ 
com UTR-ratio upon Cirbp or Rbm3 depletion, respectively 
(Supplementary Figs. S4a and S4b, Supplementary Table S8). 
These results indicated that 3 'UTRs tended to be shortened upon 
Cirbp or Rbm3 depletion. 

To examine the relationship of Cirbp- or Rbm3 -binding sites and 
the PAS-usage switching, we compared the Cirbp- or Rbm3-binding 
sites on the tandem UTRs with the ext/com UTR-ratio changes upon 
Cirbp or Rbm3 depletion (Figs. 4c and 4d). Among the APA events 
that contain Cirbp-binding sites in the common 3'UTRs (Fig. 4c), 
182 events showed downregulated ext/com UTR-ratio (Proportion 
test, P < 1.0E-15), and only 47 events showed upregulated ext/com 
UTR-ratio upon Cirbp depletion under the criteria (Fisher's exact 
test, P < 0.01; the ext/com UTR-ratio change > 0.6). If Cirbp bound 
to the extended 3'UTRs, there were few PAS-usage switching events 
upon Cirbp depletion (61 downregulated and 57 upregulated). 
Similarly, the APA events containing Rbm3 binding sites in the 
common 3'UTRs were more likely to show downregulated (78 
events; Proportion test, P < 1.0E-15) than upregulated ext/com 
UTR-ratio (8 events) upon Rbm3 depletion (Fig. 4d), whereas with 
Rbm3-binding sites in the extended 3'UTRs, there were rare PAS- 
usage switching events upon Rbm3 depletion (20 downregulated and 
9 upregulated). These results suggested that Cirbp and Rbm3 
repressed the usage of proximal PASs by binding to the common 
3'UTRs. The selected genes which showed PAS-usage switching 
upon Cirbp or Rbm3 depletion were validated by qPCR (Fig. 4e 
and Supplementary Fig. S4g). 
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Figure 4 | RNA-seq and PAR-CLIP analyses indicated that Cirbp andRbm3 repressed the usage of proximal PASs by binding to the common 3'UTRs. 

(a) Scatterplot of the common and extended 3'UTR expression of tandem UTR genes upon Cirbp depletion. Significantly increased expression of 
isoforms resulting from usage of distal (blue) and proximal PAS (red) are highlighted. Inset: total number of the different changes (Binomial test, 
P = 3.4E-08). (b) Scatterplot of the common and extended 3'UTRs expression of tandem UTR genes upon Rbm3 depletion as in (a) (Binomial test, 
P < 1.0E-15). (c) Histogram showed the ext/com UTR-ratio changes of the genes with Cirbp-binding sites in the tandem UTRs upon Cirbp depletion 
(Top, common 3'UTRs bound by Cirbp; Bottom, extended 3'UTRs bound by Cirbp). The PAS-usage switching events with up (blue) or downregulated 
(red) ext/com UTR-ratio were highlighted. Inset: the total number of cases in each region, (d) Histogram showed the ext/com UTR-ratio changes of the 
genes with Rbm3-binding sites in the tandem UTRs upon Rbm3 depletion as in (c). (e) QPCR analyses of the selected genes that showed downregulated 
ext/com UTR-ratio upon Cirbp or Rbm3 depletion. The statistical significance was calculated against the mock cells. See also Fig. S4G. Data are 
represented as mean +/- SEM (n = 3); Student's f-test * P < 0.05, ** P < 0.01. 



Cirbp and Rbm3 are important in the APA regulation upon cold 
shock and in the circadian clock. Considering that the expression of 
Cirbp and Rbm3 were significantly upregulated upon cold shock 
(Supplementary Fig. S5a), we examined the PAS-usage switching 
events in response to cold shock. Under the statistical criteria 
(Fisher's exact test, P < 0.01; the ext/com UTR-ratio change > 
0.6), 325 genes showed significant PAS-usage switching upon cold 
shock, and 63.4% (206) of them showed upregulated ext/com UTR- 
ratio, indicating that 3'UTRs tended to be lengthened upon cold 
shock (Fig. 5a). Considering that Cirbp and Rbm3 were upregu- 
lated upon cold shock, we supposed that the genes regulated by 
Cirbp or Rbm3 should show opposite direction of PAS-usage 
switching upon Cirbp or Rbm3 depletion comparing to cold shock. 
We identified the genes which showed opposite direction of 



PAS-usage switching between upon cold shock and RBP depletion 
(Cirbp: 129 genes; Rbm3: 123 genes; the ext/com UTR-ratio change 
> 0.4). Among these genes, 98 genes for Cirbp and 110 genes for 
Rbm3 showed 3'UTR lengthening upon cold shock and shortening 
upon RBP depletion (Fig. 5b). The selected genes that exhibited 
opposite PAS-usage switching were validated by qPCR (Supple- 
mentary Fig. S5b). 

To examine the relationship of Cirbp- and Rbm3-binding sites 
and the PAS-usage switching upon cold shock, we compared the 
Cirbp- or Rbm3-binding sites on the tandem UTRs with the ext/ 
com UTR-ratio changes upon cold shock (Figs. 5c and 5d). Among 
the APA events that contain Cirbp-binding sites in the common 
3'UTRs (Fig. 5c), 174 events showed upregulated ext/com UTR- 
ratio (Proportion test, P < 1.0E-15), and only 91 events showed 
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downregulated ext/com UTR-ratio upon cold shock under the cri- 
teria (Fisher's exact test, P < 0.01; the ext/com UTR-ratio change > 
0.6). In comparison, if Cirbp bound to the extended 3'UTRs, there 
were few PAS-usage switching events upon cold shock (67 upregu- 
lated and 28 downregulated). Similarly, the APA events containing 
Rbm3 -binding sites in the common 3'UTRs were more likely to 
show upregulated ext/com UTR-ratio (125 events; Proportion test, 
P < 1.0E-15) than downregulated ext/com UTR-ratio (35 events) 



upon cold shock, whereas with Rbm3 -binding sites in the extended 
3'UTRs, there were rare PAS-usage switching events upon cold 
shock (25 upregulated and 5 downregulated, Fig. 5d). These results 
suggested that 3'UTRs tend to be lengthened upon cold shock and 
upregulated Cirbp and Rbm3 play an important role in this process 
by binding to the common 3'UTRs and repressing the usage of 
proximal PASs. Several examples of depth-of-coverage profiles in 
RNA-seq were shown (Supplementary Figs. S5c-S5e). 
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Figure 5 | Cirbp and Rbm3 are important in the APA regulation upon cold shock and in the circadian clock, (a) Scatterplot of the common and 
extended 3'UTRs expression of tandem UTR genes between the mock at 32°C (Cold) and at 37°C (Normal). Significantly increased expression of 
isoforms resulting from usage of distal (blue) and proximal PAS (red) are highlighted. Inset: total number of the different changes (Binomial test, 
P = 8.0E-07). (b) Heatmap of the genes that showed opposite direction of PAS-usage switching between upon cold shock and RBP depletion. 

(c) Histogram showed the ext/com UTR-ratio changes of the genes with Cirbp binding sites in the tandem UTRs upon cold shock. The PAS-usage 
switching events with up (blue) or downregulated (red) ext/com UTR-ratio were highlighted. Inset: the total number of cases in each region. 

(d) Histogram showed the ext/com UTR-ratio changes of the genes with Rbm3 binding sites in the tandem UTRs upon cold shock as in (c) . (e) Scatterplot 
of the common and extended 3 'UTR expression of tandem UTR genes between CT6 and CT18 in circadian mouse liver as in (a) (Binomial test, P = 0.01). 
(f) Heatmap of the circadian oscillating ext/com UTR-ratio of 1153 genes in mouse liver. The genes that contained Cirbp- or Rbm3-binding sites in the 
tandem UTRs were indicated. See also Fig. S6A. (g) QPCR analyses of the ext/com UTR-ratio of the selected genes in circadian mouse liver. Data are 
represented as mean +/— SEM (n = 3). 
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Because the mouse body temperature oscillates daily with a trough 
at the CT6 phase and a peak at the CT18 phase, we quantified the 
expression changes in the common and extended 3'UTRs between 
these two time points based on the circadian microarray data of 
mouse liver (Gene Expression Omnibus, GSE11923) 45 . Out of the 
72 genes that showed significant PAS-usage switching (Fisher's exact 
test, P < 0.01; the ext/com UTR-ratio change > 0.6), 63.9% of them 
showed upregulated ext/com UTR-ratio (indicating 3'UTR length- 
ening) at the CT6 phase (Fig. 5e). Furthermore, to estimate the PAS- 
usage in a circadian manner, we calculated the ext/com UTR-ratio 
for each gene in the microarray every one hour for two days 45 , and 
identified 1153 genes that exhibited circadian oscillations of the ext/ 
com UTR-ratio (COSOPT, P < 0.01, Fig. 5f), including many core 
circadian genes, such as Clock, Rorc and Tef (Supplementary Table 
S9). Moreover, we observed that the circadian peak times of these 
genes were significantly enriched in the CT6 phase (236, 40.1%; 
Proportion test, P < 1.0E-15; Supplementary Fig. S6a). Interest- 
ingly, we found that a significant proportion of these genes contain- 
ing Cirbp- or Rbm3-binding sites in the tandem UTRs (164 and 104; 
Fisher's exact test, P = 5.0E-4 and 1.4E-5, respectively). These results 
indicated that the usage of PASs was under strong circadian regu- 
lation, and longer 3'UTRs were preferred at low body temperature. 
Cirbp and Rbm3 are potential regulators in this process. In addition, 
we estimated the PAS-usage of the genes which were regulated by 
Cirbp or Rbm3 (Fisher's exact test, P < 0.01; the ext/com UTR-ratio 
change > 0.6) in circadian mouse liver, and found that 35 and 23 
genes regulated by Cirbp and Rbm3, respectively, showed strong 
circadian oscillations of the ext/com UTR-ratio (COSOPT, P < 
0.01, Supplementary Figs. S6b and S6c). The selected genes that 
exhibited circadian oscillations of the ext/com UTR-ratio in mouse 
liver were validated by qPCR (Fig. 5g). We verified that two genes, 
Tardbp (Tdp43) and Maprel, which contained Cirbp- or Rbm3- 
binding sites in the common 3'UTRs, showed 3'UTR lengthening 
upon cold shock and 3'UTR shortening upon Cirbp or Rbm3 deple- 
tion, and exhibited strong circadian oscillations of the ext/com UTR- 
ratio in mouse liver (Figs. 6a and 6b). Collectively, these results 
suggest a novel posttranscriptional mechanism that Cirbp and 
Rbm3 regulate the amplitude of temperature entrained circadian 
gene expression by controlling APA. 

Discussion 

Temperature has long been suggested as an important entrainment 
cue in the circadian clock through posttranscriptional regulation 46 . 
For example, in Neurospora, temperature changes induce a ther- 
mosensitive alternative splicing of a core circadian gene FRQ 
(FREQUENCY) 47 . In Drosophila, a thermosensitive splicing event 
in the 3'UTR of per (period) plays an important role in temperature 
compensation and seasonal adaptation 48 . However, in mammals, the 
posttranscriptional mechanisms in the temperature entrainment of 
circadian clock remain unclear. 

In this study, we identified two cold-induced RBPs, Cirbp and 
Rbm3, as important regulators in the temperature entrained cir- 
cadian gene expression, and proposed a novel posttranscriptional 
mechanism in this process. The depletion of Cirbp or Rbm3 signifi- 
cantly reduced the oscillating amplitudes of core circadian genes in 
the temperature synchronized cells. With a combination of PAR- 
CLIP and RNA-seq analyses, we observed that the 3'UTR binding 
sites of Cirbp and Rbm3 were enriched within 100 nucleotides 
upstream of the polyadenylation sites, and the two RBPs repressed 
the usage of proximal PASs by binding to the common 3'UTRs. 
Furthermore, in mouse liver, we estimated the PAS-usage of the 
genes which were regulated by Cirbp or Rbm3, and found that many 
of them indeed showed strong circadian oscillations. In summary, 
our results indicate that Cirbp and Rbm3 modulate the amplitude of 
temperature entrained circadian gene expression by controlling APA. 



Based on these results, we proposed a model for the posttranscrip- 
tional mechanism in the temperature entrainment of circadian clock 
(Fig. 6c). In mammals, the master clock in SCN sets the daily body 
temperature cycle that affects the peripheral clocks. In the low- 
temperature phase, the cold-induced RBPs (Cirbp and Rbm3) are 
upregulated. They repress the usage of the proximal PASs of their 
target genes through binding to the common 3'UTRs, thus lengthen 
the 3'UTR and increase the posttranscriptional cis elements 28 . In 
the high-temperature phase, the downregulation of cold-induced 
RBPs leads to the preferred use of the proximal PASs and the short- 
ening of 3'UTR, thus increases mRNA stability through reducing 
microRNA-mediated repression 23 . The APA regulation in the cir- 
cadian clock could control the mRNA stability, translation and local- 
ization of the circadian oscillating genes, and influence their 
circadian functions. 

RBPs were shown to cross-regulate each other 49 . To assess the 
RBPs which are the targets of Cirbp and Rbm3, we obtained 415 
mouse RBP genes from the RBPDB database 50 , and found that 199 
of them (Fisher's exact test, P < 1.0E-15) were bound by either Cirbp 
or Rbm3. This result is consistent with the GO category enrichment 
of RNA binding in the Cirbp- or Rbm3-target genes. Moreover, 
RNA-seq analyses showed that 56 Cirbp- or Rbm3-target RBP genes 
(Fisher's exact test, P < 1.0E-15) exhibited expression changes at 
either transcript (log2 [fold change] > 0.8) or APA level (Fisher's 
exact test, P < 0.01; the ext/com UTR-ratio change > 0.6) upon 
Cirbp or Rbm3 depletion (Supplementary Fig. S6d). Intriguingly, 
several Cirbp- or Rbm3 -target RBP genes were responsible for poly- 
adenylation and 3'UTR processing, including Cpsf6, Cpsf7, 
Hnrnphl, Pabpnl and Ptbpl. In particular, Hnrnphl was shown to 
bind to Cirbp and Rbm3 mRNA, and stabilized their mRNA level 49 . 
We found Hnrnphl was bound by Cirbp and upregulated upon 
Cirbp depletion, indicating that Hnrnphl and Cirbp/Rbm3 formed 
a regulatory feedback loop. Moreover, we found that Rbm3 was 
bound by Cirbp, and the expression of Rbm3 was upregulated upon 
Cirbp depletion. Analysis of the mRNA sequences of Cirbp and 
Rbm3 suggests that while both of them are highly conserved across 
a wide range of species, Rbm3 is only present in mammals and is 
likely a result of gene duplication of Cirbp in mammals (data not 
shown). Future work should reveal the regulatory and evolutionary 
relationship of these two RBPs. In addition, another Cirbp-target 
gene, Sfpq (Psf), was shown to play a critical role in regulating cir- 
cadian clock through the recruitment of Sin3a and rhythmically 
delivering histone deacetylases to the Perl promoter to repress its 
transcription 51 . Therefore, these findings suggest that Cirbp and 
Rbm3 reside in an interconnected RBP network that regulates the 
RNA metabolism and the circadian clock. 

For the Cirbp- or Rbm3-binding sites, in addition to falling in the 
tandem UTRs, we also identified more than 3000 binding sites 
located in the CDS region and non-tandem UTRs. Besides, our 
RNA-seq analyses revealed that Cirbp and Rbm3 stabilized their 
target transcripts. These results implied that Cirbp and Rbm3 may 
have other mechanisms in regulating the circadian clock besides the 
APA regulation, such as modulating the mRNA localization and sta- 
bility of their associated circadian genes. Morf et al. recently reported 
that Cirbp is required for high- amplitude circadian gene expression. 
Using a biotin-streptavidin-based cross-linking and immunoprecipi- 
tation (CLIP) procedure, they found that Cirbp-binding sites are 
significantly enriched near the PASs. They found that Cirbp bound 
several transcripts encoding core circadian oscillator proteins, includ- 
ing CLOCK. They surmised that Cirbp influenced the circadian clock 
through increasing the cytoplasmic accumulation of mRNAs encod- 
ing Clock and other circadian regulators 21 . Collectively, these results 
suggested that Cirbp and Rbm3 regulated circadaian clock in diverse 
posttranscriptional mechanisms. 

For the core circadian genes whose ext/com UTR-ratios showed 
strong circadian oscillations in mouse liver, their ext/com UTR-ratios 
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seem not influenced by the depletion of Cirbp or Rbm3 in the cells. 
This is perhaps due to the differences between tissue and cell line, or 
there are still other RBPs which may have potential functions in the 
APA regulation of circadian clock. As shown in Table 1, several other 
RBPs such as Fus and Hnrnpdl, also showed oscillating expression 
patterns in more than six different mouse tissues, and independent 
expresion of the local molecular clock. Further studies on the func- 
tion and interaction of these remaining RBPs and the detailed 
mechanisms of APA regulation on the circadian genes will provide 
valuable information on the posttranscriptional regulation in the 
mammalian circadian clock. 

Methods 

Animals. Male BALB/c mice were housed in 12 h light/12 h dark (LD) cycles for two 
weeks, then subjected to complete darkness (DD) as a continuation of the dark phase 
of the last LD cycle. Mice were sacrificed from CT18 every four hours during two days 



in continuous dark phase. Tissues were immediately frozen in liquid nitrogen and 
stored at — 80'C until the extraction of total RNA. The age of the animals was between 
3 and 4 months. The use and care of animals and the experiments complied with the 
guidelines of the Animal Advisory Committee at the Shanghai Institutes for 
Biological Sciences, Chinese Academy of Sciences. 

Cell culture. Immortalized MEFs and 293 T cells were cultivated at 37 L 'C with 5% 
C0 2 in high glucose Dulbecco's modified Eagle's medium (DMEM) with 10% fetal 
bovine serum (FBS) (Hyclone) and antibiotic (100 U/ml penicillin and 100 ug/ml 
streptomycin). Cell lines stably expressing Flag-tagged proteins were generated by 
retroviruses co-transfection. For the treatment of cold shock, the cells were cultivated 
at 32°C for 12 h. Prior to harvest, cells were rinsed with ice-cold PBS twice. 

Square wave temperature synchronization model. 5 X 10 5 MEFs were plated into 
six cm dishes four days before the treatment was applied. Once the MEFs were fully 
confluent (90-95%), the temperature cycles that are 12 h at 37°C and 12 h at 33°C 
were used to entrain the circadian clock of the cells. After six days of entrainment, we 
either kept the temperature at constant 37"C, or continued the temperature cycles. 
The cells were harvested every four hours from the beginning of the sixth 33 "C during 
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two days. At indicated time, the cells were washed twice with ice-cold PBS and 
harvested in 1 ml of TRIzol reagent (Invitrogen). These samples were frozen and 
stored at -80°C until the extraction of total RNA. 

Accession number. The Gene Expression Omnibus accession number for the PAR- 
CLIP and RNA-seq raw data reported in this paper is GSE40468. 

Retrovirus plasmids and infection, Western blot and antibodies, RNA extraction 
and qPCR, Tandem UTRs database construction, RNA-seq, RIP and qPCR, 
Microarray analyses, PAR-CLIP library construction, sequencing and analyses are 
available in the Supplementary information. 
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